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Preface 


For  many  applications  involving  laser  irradiance  on  a  target,  there  exists  a  limitation  on 
intensity  because  of  atmospheric  effects.  One  of  the  main  atmospheric  influences  is  the 
effect  due  to  the  change  in  index  of  refraction  caused  by  temperature  gradients  on  the 
order  of  0.1  to  1°  C.  In  2007,  the  U.S.  Army  Research  Laboratory  (ARL) 
Survivability/Lethality  Analysis  Directorate  (SLAD),  White  Sands  Missile  Range 
(WSMR),  NM  decided  to  investigate  this  effect  with  Dr.  Voelz,  a  faculty  member  at  New 
Mexico  State  University  (NMSU),  Las  Cruces,  NM.  Dr.  Voelz  has  extensive  experience 
with  atmospheric  turbulence  effects  and  the  propagation  of  light  through  turbulent  media. 
Other  effects  that  are  caused  by  propagation  through  the  atmosphere  include  absorption, 
scattering,  and  high  power  laser  effects  such  as  blooming.  However,  the  effect  of 
turbulence  was  the  main  consideration  in  this  report.  The  NMSU  group  applied  both 
analytic  expressions  and  propagation  modeling  tools  in  order  to  derive  the  intensity 
profile  at  the  target  area.  Various  propagation  distances  as  well  as  initial  beam  exit 
apertures  were  considered.  The  wavelength  range  and  the  propagation  angle  were  also 
varied.  Specific  parameters  related  to  location  and  turbulence  strength  were  supplied  as 
needed.  This  effort  fills  a  gap  in  the  existing  data  already  available  for  the  propagation  of 
laser  light  through  turbulence  and  the  resulting  intensity  profiles.  Dr.  Joseph  Montoya,  of 
ARL/SLAD,  created  the  statement  of  work  and  gave  guidance  at  various  times  during  the 
project. 
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1.  Introduction 


For  many  applications  involving  laser  propagation  through  the  atmosphere,  there  exists  a 
limitation  on  the  received  flux  because  of  atmospheric  effects.  One  of  the  main  atmospheric 
influences  is  the  effect  due  to  the  change  in  index  of  refraction  caused  by  temperature  gradients 
on  the  order  of  0.1  to  1°  C.  These  gradients  comprise  the  effect  known  as  turbulence,  which 
increases  the  beam  divergence  and  creates  a  temporal  variation  of  the  beam  irradiance  at  a 
receiver.  This  variation  is  commonly  termed  “scintillation”.  For  many  communication  and 
remote  sensing  applications,  a  beam  needs  to  be  transmitted  near  the  Earth’s  surface  along  a 
path  that  is  horizontal  or  has  a  slight  incline  to  the  surface.  This  inclined  path,  or  “slant”  path, 
involves  the  transit  of  the  beam  through  varying  turbulence  strengths  at  different  altitudes  above 
the  surface. 

In  this  report,  we  list  the  currently  favored  analytical  expressions  that  describe  the  irradiance 
profile  and  scintillation  index  for  a  Gaussian  beam  propagating  through  turbulence.  These 
expressions  have  been  applied  extensively  to  free-space  laser  communications  and  other  beam 
pointing  applications.  The  turbulence  strength  Cn2  model  (daytime  and  nighttime  versions),  is 
also  discussed  and  is  generally  consistent  with  our  understanding  of  the  conditions  at  White 
Sands  Missile  Range  (WSMR),  NM.  The  expressions  and  turbulence  strength  model  are  used  to 
generate  a  set  of  plots  and  curves  that  show  the  received  beam  size  and  scintillation  index  as  a 
function  of  the  initial  beam  size  (radius)  wo,  the  wavelength  X,  the  path  elevation  angle  8,  and 
the  propagation  distance.  The  range  of  values  considered  are  as  follows:  wq  (2.5  to  10  in.), 

X  (0.5  to  5  pm),  6  (0  to  40°),  and  distance,  or  range  (0  to  20  km).  The  curves  can  be  used  to 
determine  the  expected  beam  size  and  scintillation  index  for  various  horizontal  and  slant  path 
situations. 


2.  Analytical  Expressions 


Figure  1  shows  a  schematic  diagram  for  a  laser  beam  propagating  through  the  atmosphere  to  a 
target  along  (a)  a  horizontal  path,  and  (b)  a  slant  path.  Consider  a  Gaussian  beam  with  the 
optical  field  of  unit  amplitude  written  as 
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where  wo  defines  the  initial  beam  size  (e'1  field  radius),  x  and  y  are  the  horizontal  and  vertical 
coordinates  of  the  incident  beam  field  from  the  beam  center,  respectively.  As  the  Gaussian  beam 
propagates  over  a  distance  L,  the  characteristics  of  the  received  beam  have  been  carefully 
described  by  Andrews  and  Phillips  ( 1 ,  2).  The  relevant  expressions  for  the  horizontal  and  slant 
propagations  are  now  reviewed. 


Figure  1.  Schematic  diagram  of  a  Gaussian  beam  propagating  through  turbulence 
along  a  (a)  horizontal  and  (b)  slant  path. 

2.1  Horizontal  Propagation 


2.1.1  Intensity  Profile 

The  irradiance  (often  referred  to  as  intensity)  function  of  the  beam  at  a  distance  L  from  the 
source  can  be  expressed  by 


wn 
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where  w{L)  is  the  average  beam  size  (radius)  at  the  receiver  given  by 
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where  k  =  2nl 'X  and  X  is  the  wavelength  of  the  beam,  or  is  the  Rytov  variance,  and  A(L)  is  the 
Fresnel  ratio  for  vacuum  propagation. 

The  Rytov  variance  ctr  is  the  variance  of  the  log-irradiance  fluctuations  of  the  field.  It  generally 
signifies  “weak”  turbulence  (phase  dominated  effects)  if  it  has  a  value  less  than  1,  and  “medium- 
to-strong”  turbulence  if  greater  than  1.  If  the  turbulence  strength  along  a  horizontal  path  is 
assumed  to  be  constant  then  the  Rytov  variance  can  be  expressed  as 

7  n 

<r2R(L)  =  1.23C2Hk*lJ,  (4) 

where  C„  is  the  turbulence  structure  constant  that  characterizes  the  strength  of  the  index  of 
refraction  fluctuations.  Equation  4  assumes  a  Kolomogorov  turbulence  spectrum.  Note  that  the 
Rytov  variance  is  directly  proportional  to  the  turbulence  strength,  approximately  inversely 
proportional  to  the  wavelength  and  approximately  proportional  to  the  propagation  distance 
squared. 


The  Fresnel  ratio  A(L)  in  equation  3,  characterizes  the  beam  spread  due  to  diffraction  for  vacuum 
propagation.  A (L)  is  associated  with  the  beam  curvature  parameter  0(Z),  which  describes  the 
phase  curvature  of  the  beam  as  it  propagates  in  vacuum.  These  parameters  are  dependent  on  the 
initial  Fresnel  ratio  A0(T)  and  curvature  parameter  0O (L)  that  include  the  initial  beam  size  and 
curvature  at  the  transmit  plane.  The  key  relationships  are  as  follows: 
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where  F0  is  the  (phase)  curvature  of  the  beam  at  the  transmitter.  For  a  Gaussian  beam  F{)  is 
usually  assumed  to  be  infinite,  so  ©o (L)  =  1. 

2.1.2  Scintillation  Index 

The  scintillation  index  is  defined  as  the  normalized  variance  of  the  irradiance  fluctuations  at  a 
receiver,  that  is, 

ll\x,y,L))-(l{x,y,L))2 

aj(x,y,L)  = - - - -2 - ,  (9) 

(l(x,y,L))2 

where  / is  the  irradiance,  x  and}’  are  transverse  coordinates  at  the  receiver,  and  L  is  again  the 
propagation  distance.  Scintillation  index  expressions  are  commonly  separated  into  two  forms: 
one  for  “weak”  fluctuation  conditions  and  one  for  “strong”  fluctuation  conditions.  As  previously 
noted,  the  test  for  these  conditions  involves  the  Rytov  variance. 


Based  on  a  Kolmogorov  turbulence  spectrum,  the  scintillation  index  at  the  beam’s  center  point 
for  a  weak  turbulence  case  can  be  expressed  as 


iw(x,y,L)  =  3.86 o\  Re  i«f(®(L)  +  iA(Z))-^(A(Z))6 
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where  F(. . .)  is  a  hypergeometric  function  (1),  ©  (Z)is  equal  to  1-  Q(L)  and  i  is  the  imaginary 
unit.  For  a  Gaussian  beam,  the  initial  beam  curvature  Fo  is  oo,  and  equation  10  becomes 

^  =  3.86a2  jo.4[(l  +  20(Z))2  +4(A(T))2]^  cos^tan'1  f 1  ± ^®(Z,)  jl - ^(A(Z))f  } ■  (11) 
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Under  medium-to- strong  turbulence  conditions,  the  scintillation  index  a \  s  is  given  by 

0.49cr7  w2  0.5  la^2 

7  5 

12  Afi  (  12  he 
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J  \  J 

where  the  scintillation  index  for  the  weak  case  a/  ,,  2  is  applied  in  equation  12  to  determine  the 
result  for  the  medium-to-strong  case. 
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2.2  Slant  Path 


2.2.1  Intensity  Profile 

For  Gaussian  beam  propagation  through  turbulence  over  a  slant  path  of  length  L,  the  expression 
of  the  received  beam  size  is  generally  the  same  as  equation  3.  However,  in  this  situation,  the 
path  is  characterized  by  the  elevation  angle  #(i.e.,  the  angle  of  beam  with  respect  to  the 
horizontal  direction),  and  the  height  change  //between  the  transmitter  and  receiver.  For  an 
upward  slant  path,  the  relationship  between  these  parameters  can  be  easily  obtained  from  the 
geometry  in  figure  1(b): 


H  =  L  sin  (9.  (13) 

Because  the  beam  height  changes  along  the  propagation  path  and  the  turbulence  strength  Cn2  is  a 
function  of  height,  the  Rytov  variance  is  modified  to  the  following  form 


=  2.25k 6  H  6  (sin  0)~T  ^+h°  C2  (h) 


h) 


h-h0 

H 


dh , 


(14) 


where  C„2(h)  indicates  the  turbulence  strength  as  a  function  of  height  h  from  the  ground  and  ho  is 
the  initial  height  of  the  beam  transmitter.  Substituting  equation  14  into  equation  3  gives  an 
expression  for  the  beam  size  for  the  upward  slant  path  propagation. 


2.2.2  Scintillation  Index 

The  expressions  within  the  curly  bracket  of  equation  1 1  involve  only  the  parameters  0  and  A, 
which  are  functions  of  propagation  distance  L,  but  are  independent  of  the  turbulence  parameters. 
Consequently,  the  Rytov  variance  ctr  is  the  only  parameter  affected  by  the  slant  path  turbulence 
(2).  Using  equation  14  with  equation  11,  leads  to  an  expression  for  the  scintillation  index  for 
slant  path  propagation  under  weak  turbulence.  The  scintillation  index  for  medium-to-strong 
turbulence  is  again  given  by  equation  12. 


3 .  C2  Profiles 


The  C2  profile  is  the  key  function  that  characterizes  the  turbulence  strength  with  height. 
Information  provided  by  WSMR  personnel  suggested  that  typical  daytime  conditions  at  WSMR 
can  be  modeled  with  a  -4/3  height  dependence  near  the  ground  and  nighttime  conditions  involve 
a  -2/3  dependence.  A  C2  profile  published  by  Comeron  et  al.,  matches  these  conditions; 
therefore,  we  used  this  profile  for  the  results  presented  in  this  report  ( 3 ). 
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The  expression  used  for  the  height  dependence  of  C„  for  daytime  conditions  is 

4  2 
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where  the  surface  layer  height  hs  is  4  m  and  the  reference  height  ho  is  8  m.  The  height  marking 
the  fall-off  of  the  h 4/3 dependence  h,  is  100  m,  the  height  above  the  terrain  is  h,  and  the  height 
characterizing  the  exponential  fall-off  of  the  structure  constant  in  the  free-atmosphere  region  h\  is 
1500  m.  The  tropopause  height  ht  is  13,000  m;  the  root  mean  squared  wind  velocity  averaged 
over  the  5  to  20  km  altitude  interval  w\  is  30  m/s;  the  structure  constant  in  the  zone  with  h~ 4/3 
dependence  is  C„o  =1.7x10"  m"  ;  and  the  structure  constant  value  for  free  atmosphere  is  Cn\ 

=  10"16  m"2/3, 

The  expression  used  for  the  height  dependence  of  Cn2  for  nighttime  conditions  is 
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if  h>  ht 


where  h,  =30  m,  Cno2  =  9.5x  10"15m"2/3,  and  Cn2  =  4.5  x  1 0‘l7m‘23.  Figure  2  presents  a  plot  of  the 
logarithm  of  the  C„  profile  for  day-  and  nighttime  conditions.  It  is  clear  that  from  the  ground  up 
to  about  100  m  the  daytime  turbulence  is  significantly  stronger  than  nighttime  turbulence. 
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Figure  2.  Daytime  and  nighttime  C„2  profiles  used  in  this  report  (5). 


4.  Results  and  Discussions 


4.1  Beam  Size 

For  our  calculations,  the  initial  height  of  the  laser  transmitter  ho  was  assumed  to  be  8  m.  A  set  of 
curves  for  the  beam  size  as  a  function  of  propagation  distance  L,  elevation  angle  6,  initial 
Gaussian  beam  size  w0,  wavelength  X,  and  different  C,,2  profiles  are  shown  in  figures  3  through 
6.  Observations  regarding  these  results  are  as  follows: 

1.  For  horizontal  propagation,  the  beam  size  is  driven  by  the  atmospherics  so  wavelength  and 
initial  beam  size  have  little  effect  on  the  beam  divergence  where  the  turbulence  is  most 
severe. 

2.  Diffraction  limited  spreading  becomes  more  of  a  factor  for  the  slant  path  cases  where  the 
turbulence  is  less  severe,  thus  the  wavelength  and  initial  beam  size  become  more 
important. 

3.  Beam  spread  is  somewhat  larger  during  the  day  than  at  night,  but  the  change  in  size 
between  these  two  cases  is  not  substantial. 
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Figure  3.  Beam  size  w(Z)  for  the  slant  path  propagation  of  a  Gaussian  beam  versus  propagation  distance  L  in 
the  daytime. 


Note:  X  =  0.5  pm  in  (a),  (b),  and  (c),  1  pm  in  (d),  (e),  and  (f),  elevation  angle  6  =  0°,  10,  40°,  and  initial 
Gaussian  beam  size  w0  =  2.5  in.  (solid),  5  in.  (dot),  7.5  in.  (dash),  and  10  in.  (dash-dot). 
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Figure  4.  Beam  size  w(L)  for  the  slant  path  propagation  of  a  Gaussian  beam  versus  propagation 
distance  L  in  the  daytime. 

Note:  with  X  =  2  pm  in  (a),  (b),  and  (c),  5  pm  in  (d),  (e),  and  (f),  elevation  angle  6  =  0°,  10,  40°,  and  initial 
Gaussian  beam  size  w0  =  2.5  in.  (solid),  5  in.  (dot),  7.5  in.  (dash),  and  10  in.  (dash-dot). 
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Figure  5.  Same  as  figure  3,  except  in  the  nighttime. 
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Figure  6.  Same  as  figure  4,  except  in  the  nighttime. 
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4.2  Scintillation  Index 


A  set  of  curves  for  the  scintillation  index  are  presented  in  figures  7  through  10.  The  observations 
on  these  results  are  as  follows: 

1.  As  expected,  longer  wavelengths  and  higher  elevation  angles  generally  result  in  lower 
scintillation. 

2.  For  horizontal  propagation,  the  scintillation  index  saturates  (at  a  value  of  ~1.7)  for  all 
wavelengths  and  initial  beam  sizes  within  the  20  km  distance,  although  the  shortest 
wavelengths  reach  saturation  with  propagation  lengths  of  less  than  1  km. 

3.  For  the  slant  paths,  the  smaller  initial  beam  sizes  result  in  somewhat  less  scintillation. 
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Figure  7.  Scintillation  index  versus  propagation  distance  L  during  daytime. 

Note:  X  =  0.5  pm  in  (a),  (b),  and  (c),  1  pm  in  (d),  (e),  and  (f),  elevation  angle  6  =  0°,  10,  40°,  and  initial 
Gaussian  beam  size  w0  =  2.5  in.  (solid),  5  in.  (dot),  7.5  in.  (dash),  and  10  in.  (dash-dot). 
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Note:  X  =  2  jam  in  (a),  (b),  and  (c),  5  jam  in  (d),  (e),  and  (f),  elevation  angle  0  =  0°,  10,  40°,  and  initial 
Gaussian  beam  size  w0  =  2.5  in.  (solid),  5  in.  (dot),  7.5  in.  (dash),  and  10  in.  (dash-dot). 
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Figure  9.  The  same  as  figure  7,  except  in  the  nighttime. 
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Figure  10.  The  same  as  figure  8,  except  in  the  nighttime. 
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5.  Concluding  Remarks 


The  curves  provided  in  this  report  are  adequate  for  first-order  estimates  of  beam  size  and 
scintillation  index  for  the  parameter  space  that  was  investigated.  Some  of  the  results  were  spot- 
checked  with  results  from  wave  optics  simulations  and  were  found  in  good  agreement.  In  regard 
to  beam  size,  there  are  some  indications  from  our  wave  optics  simulation  work,  that  beam  size  in 
strong  turbulence  may  be  overestimated  by  the  analytic  theory  applied  in  this  report.  There  may 
be  a  saturation  effect  that  is  not  being  considered;  therefore,  this  issue  requires  further 
investigation.  We  emphasize  that  scintillation  index  values  are  computed  only  for  the  center  of 
the  received  beam  and  normalized  intensity  fluctuations  near  the  beam  edges  will  tend  to  be 
larger.  Furthermore,  the  effect  of  beam  wander  is  not  fully  characterized  in  the  analytic 
expressions  used  in  this  report  ( 4 ).  Beam  wander  becomes  important  in  situations  where  the 
beam  is  generally  tilted  rather  than  broken  up  by  turbulence.  Beam  wander  may  tend  to  increase 
the  scintillation  values  for  some  of  the  slant  path  cases.  We  did  not  investigate  the  temporal 
characteristics  of  the  scintillation,  which  depend  on  a  variety  of  factors  such  as  wind  speed.  This 
could  be  examined  in  the  future.  Finally,  we  emphasize  that  our  results  are  dependent  on  the 
assumed  C,,2  profile,  although  most  commonly  accepted  profiles  should  produce  similar  results. 
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